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I. INTRODUCTION 

The interplay between itinerant carriers in a partially filled band and localized magnetic 
moments is conventionally studied within the ferromagnetic Kondo lattice model (FKLM), 
wherein — JS/.er/ represents the local exchange interaction between the localized magnetic 
moment S/ and the itinerant electron spin 07. Magnon excitations in the FKLM provide im- 
portant information about this interplay, with the magnon dispersion determining the finite- 
temperature spin dynamics, transition temperature T c , and also providing information about 
the carrier- induced spin-spin interactions in the ferromagnetic state, negative-energy modes 
signalling instability due to competing antiferromagnetic (AF) interactions, and magnon 
damping highlighting the characteristic spin-charge coupling in a band ferromagnet. 

Magnons in the FKLM have been studied as function of electron (hole) density n (p) 
in the conduction (valence) band and the spin-fermion coupling J in the context of heavy 
fermion materials,- ferromagnetic metals Gd, Tb, Dy, doped EuX,- colossal magnetoresis- 
tive manganites (CMR),-^ 1 ^ ordered diluted ferromagnet s-£ relevant for the ordered double 
perovskite Sr 2 FeMo06,- and diluted magnetic semiconductors such as Gai_ x Mn x As in which 
impurity positional disorder and clustering play a very significant role on magnon excitations 
and hence on the finite-temperature spin dynamics.- 1 ^ 

Quantum corrections to the magnon spectrum at finite S, not accounted for in the 
above leading-order calculations, have also been studied using various approaches such as 
exact diagonalization,— ^ modified RKKY approach,— variational method,— 1 ^ the Holstein- 
Primakoff transformation.— 1 ^ Most of these theoretical investigations^ were carried out in 
the strong-coupling (double-exchange) limit (J/W ^> 1), where W is the bandwidth of 
the mobile electrons and J is their exchange coupling to the localized spins. Although 
providing a good description of the experimentally observed magnon damping in the fer- 
romagnetic phase of colossal magneto-resistive (CMR) manganites,—, these investigations, 
however, could not satisfactorily account for the characteristic anomalous magnon softening 
observed for zone-boundary modes. 

For ferromagnetic manganties, typical parameter values are: t = W/12 ~ 0.2 — 0.5 eV 
for the mobile (e g ) electrons and 2J ~ 2 eV for their exchange coupling to the localized 
core (t2g) spins,— so that with 2J/t ~ 4 — 10, the intermediate-coupling regime (J/W ~ 1) 
appears to be more appropriate. Theoretical studies have therefore been carried out recently 
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in this regime as well, such as the variational investigation,— where anomalous softening is 
pronounced for J/W ~ 1. Furthermore, by taking into account the Coulomb repulsion be- 
tween band electrons, which is the largest energy scale in manganites and often omitted in 
conventional FKLM investigations, several realistic features such as doping dependent asym- 
metry of the ferromagnetic phase and enhanced zone-boundary anomalous softening have 
been demonstrated, thereby highlighting the importance of correlated motion of electrons 
on spin dynamics.— ^ 

Spin dynamics in the ferromagnetic phase of CMR manganites has attracted consid- 
erable current interest.— Recent spin- wave excitation measurements have revealed several 
anomalous features in the magnon spectrum near the Brillouin-zone boundary.— These 
observations are of the crucial importance for a quantitative understanding of the carrier- 
induced spin-spin interactions and magnon damping, and have highlighted the limitations 
of various existing theoretical approaches. For example, the predictions of disorder- induced 
softening 26 of the zone-boundary modes and magnon-phonon coupling 27 as the origin of 
magnon damping are in sharp contradiction with these observations. Furthermore, the dra- 
matic difference in the sensitivity of long-wavelength and zone-boundary magnon modes on 
the mobile charge carriers has emerged as one of the most puzzling feature. Observed for 
a finite range of hole concentrations, while the spin stiffness remains almost constant, the 
softening and broadening of the zone-boundary modes show substantial enhancement with 
increasing hole concentrations.— 1 ^ 

It is therefore of interest to develop an approach wherein quantum fluctuation effects 
associated with both finite S and electron correlation can be studied on an equal footing. 
In this paper we will introduce a purely fermionic representation of the FKLM in terms of a 
multi-orbital Hubbard model, which allows conventional diagrammatic tools to be employed 
for the investigation of quantum corrections beyond the leading order. We will make use of 
the inverse-degeneracy expansion scheme, employed recently to investigate ferromagnetism 
in the Hubbard model,- 28 which systematically incorporates correlation effects in the form of 
self-energy and vertex corrections, while preserving the spin-rotational symmetry and hence 
the Goldstone mode explicitly. This purely fermionic analysis can be seamlessly extended 
to correlated band electrons, which will be reported separately. 

The organization of this paper is as follows. A purely fermionic representation for the 
FKLM is introduced in section II. The magnon propagator is discussed in section III in 
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terms of the irreducible particle-hole propagator, for which the first-order quantum cor- 
rections incorporating self-energy and vertex corrections are obtained in section IV within 
a systematic inverse-degeneracy expansion scheme. Results for the magnon self energy, 
renormalized magnon energy and damping, and electronic spectral-weight transfer due to 
electron-magnon coupling are presented in section V for one-, two-, and three-dimensional 
cases. Finally, conclusions are presented in section VI. 



The conventional FKLM model involves mobile electrons in a partially-filled band 
exchange-coupled to localized spins. In this section we introduce a purely fermionic rep- 
resentation in terms of a multi-orbital Hubbard model, which allows conventional diagram- 
matic tools to be employed to investigate quantum corrections beyond the random phase 
approximation (RPA). For concreteness, we start with a two-orbital Hubbard model: 



with a magnetic (a) orbital and a non-magnetic (f3) orbital. Here the correlated and localized 
a orbital with large Hubbard interaction U yields exchange-split spin bands with vanishing 
bandwidth and, with occupancies n a j = 1 and n a i = at half-filling, represents spin S — 1/2 
localized magnetic moments. Spin-S* magnetic moments, such as S = 5/2 in DMS systems 
and S = 3/2 in manganites, can be similarly represented by simply including multiple 
(A/" = 2S) a orbitals per site with strong Hund's coupling. The (3 band with dispersion 
represents mobile fermions such as valence-band holes in DMS or e g electrons in manganites. 
The Hund's coupling J between the two orbital spins S ia and represents the conventional 
exchange interaction between the localized spin and the mobile fermion spin, as considered 
phenomenologically in the ferromagnetic Kondo lattice model. A correlation term can also 
be included for the mobile (j3) fermions. 

In the saturated ferromagnetic state, with ordering chosen in the z direction, and spa- 
tially uniform magnetizations m a = 1 and mp = m for the localized and mobile fermions, 
corresponding to fully occupied a j band and partially occupied (5 j band, the interaction 
terms reduce to: 



II. TWO-ORBITAL HUBBARD MODEL 



e aQ'iacr a iacr + / J ek a k/3cr a k/3cr ~ U / J ^ia-^ia ~ 2 J ^ ^ Sjo-.Sj^ 
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FIG. 1: Diagrammatic representation of coupled equations for components of x + (shaded box) 
in terms of the irreducible particle-hole propagator 0(q, u) (open box). 

at the Hartree-Fock level, where ip^ represent the fermion field operators for the two orbitals 
\i = a, j3, and the corresponding exchange splittings are given by: 

2A a = 2(U(S? a ) + J(£fy)) = (U + Jm) 

2A P = 2J(S* a ) = J . (3) 

The bare antiparallel-spin particle-hole propagators for the two orbitals are given by: 
1 1 



X>) 
X?(q 



2A a + uj — ir\ U + Jm + uj — irj 

El m 
— — = -TT-. — — ( for q = °) > 



J + uj — it] 



(4) 



where e£ = ek — aAp are the exchange-split band energies, and superscripts + (— ) refer to 
particle (hole) states above (below) the Fermi energy ep. 



III. MAGNON PROPAGATOR 

The different components (fi,fi' = a, f3) of the time-ordered magnon propagator in the 
ferromagnetic state |\&o) are given by: 



dt e M*~0 X)e^( r '-^)(*o|T[5r(t)5+,(f )]*o> 
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(5) 



in terms of the spin-lowering and -raising operators = ipj^a^ /2)i(ji ll . As shown in Fig. 
1, the coupled equations for the different components of \~ + can De expressed exactly as: 



X 



(6) 



in terms of the irreducible particle-hole propagators 0. In analogy with the inverse- 
degeneracy (l/A/ - ) expansion studied recently for the band ferromagnet,— we consider a 





t 

a 

(a) 



(6) 



a 






< 




> 



at 



(c) 

FIG. 2: The first-order quantum corrections to the irreducible particle-hole propagator </>(q, cj) 



systematic expansion: 



f + ^) (2) 



(7) 



for the irreducible propagator, where x° are the bare particle-hole propagators, and correla- 
tion effects in the form of self-energy and vertex corrections are incorporated systematically 
in the quantum corrections ft 1 ' etc. so that spin-rotation symmetry and hence the Goldstone 
mode are explicitly preserved order by order. 
Solving these coupled equations, we obtain: 

0/3/3(1 - U(j) aa ) + U(j)p a (f) a p 



1 - 


U (j) aa - J 2 <f) aa <f)i3l3 - J (0/3a + <, 
4>/3a + J4>aa4>/3I3 ~ J^af. 


W) + J 2( t>a0(t>p a 
?0/3a 


1 - 


U 4>aa ~ J 2 4>aa4>f3(3 


- J (4>f3a + (, 
4*aa 




1 - 
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~ J {4>/3a + <, 


+ J 2 4>a/34>Pa 


g* 1 




Kl + J 2 ^ 





(8) 

X-J (9) 

(10) 

Vaa (11) 

represents the local magnon propagator for the localized spin (a orbital). Resulting from 
the coupled nature of Eq. (6) for the different components, the common denominator in 
Eqs. (8-10) ensures a single Goldstone mode for all components, as expected. 

The first-order quantum corrections to the irreducible propagator <p are shown diagram- 
matically in Fig. 2, and correspond to the following physical processes: a) positive cor- 
rection (to the bare particle-hole propagator) due to spin- J. spectral-weight transfer arising 



X/3/3 
X/3a 



Xaa 



where 
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from self-energy correction, b) negative correction due to particle-particle correlations, and 
c) positive correction due to spin- J, correlations induced by spin-f charge fluctuations. The 
corresponding expression are obtained as: 



Q k \ 

j2 E/£xiS A (Q,si)E 



+ uj — Q — IT) 

(IV. x ^ / I \ / 1 



qJ 2m k V e k-q- e k +W~iriJ V e k-q+Q- e k +CU-fi-2^ 




'' CT^m 1 ' (12) 

where 

Xrp+a(Q ' fi) = i-^(fi)-J^°(fi)xS(Q,n) (13) 

denotes the leading-order result for the magnon propagator in the RPA. 

Separating out the bare part from the quantum corrections in the irreducible particle- 
hole propagators in (10), and dropping explicitly second-order terms such as 0^0^, we 
obtain: 

Xaa + (q, [i0 + Jm _ J2 X (q? 7J}] _ [SQ(q) 7J) + S/3 ( qj 7J) + 2Sa/3 ( q? w )] ^ 

for the «a component. The zeroth-order first term in the denominator yields the RPA result. 
The Goldstone mode at this level is easily verified, confirming the systematic nature of the 
expansion. For q = 0, the bare fermion propagator x^(0,cj) = m/(J + cu), which yields a 
pole at uj = 0, as expected from the continuous spin-rotation symmetry. For finite q, the 
RPA magnon energy is obtained by solving the pole equation 

uj + Jm- J 2 x p(q,uj) = 0. (15) 

As J 2 x°(0,0) = Jm from (4), and the uj dependence of the bare fermion propagator x° is 
relatively weak, the RPA magnon energy is approximately obtained as: 

^ = ^ 2 (2S)[x°(0)-4(q)] (16) 
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for the spin-S* case. In the strong-coupling limit, and for dispersion = — 2t(cos k x +cos k y + 
cosk z ) on a simple cubic lattice, the bare magnon energy u;° ~ (3t/S) ^^n^cos k x (l — 7 q ) 
is of order 1/S and has the Heisenberg form. 

The first-order magnon self energies in Eq. (12) are obtained as: 

S Q (q,cu) = (U + Jm + u;) 2 0«(q,u;) 
E^(q,w) = J 2 ^(q,uj) 
E Q/3 (q, u) = J(U + Jm + w)0j2(q, w), (17) 

and we now discuss the resulting magnon renormalization. As the Goldstone mode is ob- 
tained at the RPA level itself, the magnon self-energy corrections must exactly cancel for 
q, uj = 0, as explicitly demonstrated below. The three magnon self energy terms, corre- 
sponding to the a, j3 and mixed contributions, can be combined in a simple expression: 

Stotal = S a + Efl + 2E a 



*0 



^ 2 E/^^(Q,fi)i:[( 



-T+ 



e k- q+Q - e k + uj-Q-ir] 



X 



1 - 



J 



k— q 



- el + u — irj 



(18) 



which identically vanishes for q, to — 0. Equation (18) highlights the spin-charge coupling 
in the magnon self energy, between magnons and charge fluctuations in the partially-filled 
majority-spin band. Magnon decay into intermediate magnon states accompanied with 
charge fluctuations results in magnon-energy renormalization and magnon damping, as dis- 
cussed below. 

It is straightforward to generalize to the case of localized spin-S* magnetic moments, 
corresponding to M — 2S localized a orbitals. Replacing the fermion-magnon interaction 
vertex J by JV2S, we obtain: 



£totai(q,") = j2 ( 2S )J2 J 



1 




2ni Q + u^ 
2JS 



17] 



1 



(19) 



In the strong-coupling (double-exchange) limit where <C 2JS, the magnon self energy 
simplifies to: 



S to tai(q^) = J 2 (2S)J2J2 

Q k 
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ek-q - e k + u; 
2JS 
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FIG. 3: The scaled self energy for the one-dimensional chain in the double-exchange limit, showing 
filling-dependent deviation from the Heisenberg form (dashed line), in agreement with the exact 
finite-size results of Ref. [12]. The deviation changes sign at n « 0.43. 
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(20) 



Q k "k-q+Q ~ c k "T- ^ -r 

This magnon self energy has the following characteristic features: i) it has the same energy 
order (t) as the bare magnon energy u^, ii) it is not of the form Xlk( e k-q — e k), implying 
deviation in the magnon dispersion from the Heisenberg form, in) it involves spin-charge 
coupling, resulting in magnon damping, iv) it is down by the factor 1/25 relative to the 
classical magnon energy, and hence vanishes in the 5 — > oo limit. Limiting behaviours of the 
resulting magnon dispersion — non-Heisenberg form in the J — > oo limit and Heisenberg 
form in the 5 — > oo limit — are in agreement with the exact calculations for the one- 
dimensional chain.— 

The imaginary part of the strong-coupling magnon self energy: 

2 

Q -£-+u + d° Q ) (21) 



-Im£ totaI (q,u;) = (— j $^( e k- q - e k + ^) 2 5(ej c + qH 

71 ^ ' Q k 



yields finite magnon damping and linewidth at zero temperature, arising from magnon de- 
cay into intermediate magnon states accompanied with charge (majority- spin particle-hole) 
fluctuations. The above result for magnon damping due to the spin-charge coupling is in 
exact agreement with earlier strong-coupling results obtained using the Holstein-Primakoff 
transformation and 1/5 expansion, but only after terms of third order in 1/5 are included, 
amounting to a self-energy correction in the second-order magnon self energy.—^ 
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FIG. 4: The scaled renormalized magnon energy u q .5 for the sc lattice in the double-exchange 
limit, showing enhanced self-energy correction with decreasing spin quantum number S. 

IV. RESULTS 

Fig. [3] shows the self energy, calculated from Eq. 20 with u = 0, for the one-dimensional 
chain in the double-exchange limit for S = 3/2. The deviation in the magnon self energy 
from the Heisenberg form is seen to be filling dependent, with the deviation changing sign 
at n m 0.43, in agreement with the exact finite-size results.— 

Fig. H] shows the scaled renormalized magnon dispersion u^.S for the simple cubic lattice 
in the double-exchange limit. There is significant magnon-energy reduction due to spin- 
charge coupling, which is enhanced with decreasing spin quantum number S highlighting 
the quantum 1/S effects, whereas it vanishes in the large S limit. Even in the strong- 
coupling limit, the magnon dispersion shows deviation from the Heisenberg form, which is 
more pronounced at finite J as discussed below. 

Fig. [5] shows the renormalized magnon dispersion for the simple cubic lattice in the 
intermediate-coupling regime. The anomalous momentum dependence of the magnon self 
energy results in significant anomalous softening which is robust even for J ~ W where 
the Stoner gap is substantial. Even for large J, the ratio of magnon energies at X and 
R is noticeably smaller than 1/3. Also seen in the spin- wave excitation measurements of 
ferromagnetic manganites,-— this feature indicates significant presence of second- {J 2) and 
third- neighbour (J3) antiferromagnetic spin-spin interactions. 

Fig. [6] shows the behavior of the magnon spectral function for a square lattice from 
the zone center to the zone boundary along the (1, 1) direction at band filling n=0.82 and 
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FIG. 5: The renormalized magnon dispersion for the simple cubic lattice in the intermediate- 
coupling regime, showing significant anomalous softening and deviation from the Heisenberg form. 
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FIG. 6: The magnon spectral function for the square lattice in the double-exchange limit at 
filling n=0.82 and for different spin quantum numbers S, showing substantial broadening of zone- 
boundary modes due to magnon damping, especially for low S. 



in the double-exchange limit. Magnon damping, as indicated by the linewidth, is seen to 
be sharply enhanced for zone boundary modes, especially for low spin quantum number S 
where quantum effects are significantly enhanced. 

Fig. [7] shows a comparison of intermediate-coupling and strong-coupling results for the 
renormalized magnon dispersion and magnon damping for a simple cubic lattice with n = 0.7 
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FIG. 7: Comparison of intermediate-coupling (a) and strong-coupling (b) results for the renormal- 
ized magnon dispersion u; q and magnon damping r q , for a simple cubic lattice with n = 0.7 and 
S = 3/2. 
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FIG. 8: The spin-resolved density of states for the simple cubic lattice at filling n = 0.5, showing 
the transfer of minority-spin spectral weight and band narrowing due to correlation effects. 

and S = 3/2. Here the magnon self energy was evaluated from Eq. (19) with — uj = uj®, the 
bare magnon energy. Relative to the magnon energy at R, magnon energies at X and M 
are substantially suppressed at finite J, highlighting the anomalous softening, while magnon 
damping is significantly enhanced near M. The strong momentum dependence of magnon 
damping, particularly along the M-R direction, highlights the role of charge fluctuations in 
the spin-charge coupling mechanism responsible for magnon damping. 
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FIG. 9: Intensity plot of the minority-spin spectral function along symmetry directions, show- 
ing a nearly coherent high-energy branch with substantially reduced band width, and a strongly 
incoherent low-energy branch with hot spots between T-X and T-R. 

Fig. [8] shows the spin-resolved density of states (DOS) of the band (f3) electrons on a 
simple cubic lattice at filling n = 0.5. For spin-j, the bare (HF) DOS is exact as there 
are no self-energy corrections in the saturated ferromagnetic state at T = 0. However, a 
spin- 1 particle can decay into a magnon and a spin-j particle due to the fermion-magnon 
interaction, as shown in the self-energy part of Fig. 2(a), resulting in considerable spin-j 
spectral- weight transfer above the Fermi energy ep. Within an approximate resummation 
procedure which incorporates particle-particle correlations,— the spin- J. self energy for the 
/3 orbital was calculated from: 



S/3|(q,w) 



y(0) 



l-Eg>(q, W )/Eg" 



(22) 



where E{jJ = 2A^ = 2 JS is the HF-level self energy, and the first-order self energy incorpo- 
rating electron-magnon interaction is given by: 

Q 



Eg(q,a;) = J 2 (2S)]T J ^ X^pa(Q^) 



u-n-el + _ Q + ir] / 
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= J 2 (2S)J2 • (23) 

q w - figj - e[_ Q + it; 

Fig. O depicts the energy-momentum dispersion of spin- j electrons for the above case in 
terms of an intensity plot (on a log scale) of the renormalized spectral function A^(u). The 
low-energy branch of the spectral-weight is strongly incoherent, with (relative) hot spots 
near the T point between T-X and T-R. The high-energy branch remains nearly coherent, 
and the dispersion retains the bare character, although the renormalized band width is 
substantially reduced. 



V. CONCLUSIONS 



In conclusion, a purely fermionic representation for the FKLM was introduced in terms 
of a multi-orbital Hubbard model involving localized and band electrons representing lo- 
cal moments and itinerant fermions, and correlations effects were investigated using a sys- 
tematic inverse-degeneracy diagrammatic scheme which explicitly preserves the continuous 
spin-rotation symmetry and hence the Goldstone mode. First-order quantum corrections to 
magnon excitations due to spin-charge coupling were investigated in the full range of inter- 
action strength, and shown to result in strong magnon damping and anomalous softening 
for zone boundary modes, which accounts for several zone-boundary anomalies observed in 
recent spin-wave measurements of ferromagnetic manganites. 

We find that in the strong-coupling limit, the magnon self energy for the one-dimensional 
chain exhibits filling-dependent deviation from the Heisenberg form, with the deviation 
changing sign near n ~ 0.43, in agreement with the exact finite-size results. 1 - For the sim- 
ple cubic lattice also, the renormalized magnon dispersion shows noticeable deviation from 
the Heisenberg form. Magnon decay into an intermediate-state magnon and majority-spin 
charge fluctuations due to spin-charge coupling lead to strong damping and line-broadening 
for zone-boundary modes in two and three dimensions, especially for low S. 

For finite J, spin-charge coupling is significantly enhanced, resulting in strong zone- 
boundary magnon damping and anomalous magnon softening even for large Stoner gap. 
This is in contrast to the RPA result^ where both these features were found only in the weak- 
coupling regime where the Stoner gap is vanishingly small and the magnon branch merges 
with the Stoner continuun. However, in this weak-coupling regime, quantum corrections are 
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strong enough to destabilize the ferromagnetic state. We also find that the ratio of magnon 
energies at X and R is noticeably smaller than 1/3, indicating significant presence of second 
( J2) an d third-neighbour (J3) antiferromagnetic spin-spin interactions. This feature is also 
seen in a recent spin-wave excitation measurement of ferromagnetic manganites. 24 

The interplay of electronic and magnetic excitations is highlighted by the fermion-magnon 
scattering of minority-spin electrons, resulting in substantial spectral-weight transfer to lower 
energies, with the additional spin-j particle - spin-f hole process being responsible for the 
quantum exchange correction to spin stiffness. Intensity plot of the renormalized minority- 
spin spectral function showed a strongly incoherent low-energy branch, with (relative) hot 
spots near the T point between T-X and T-R. The high-energy branch remained nearly 
coherent and the dispersion retained the bare character, although the renormalized band 
width is substantially reduced. 

Finally, this purely fermionic representation and diagrammatic analysis not only provides 
insight into quantum effects on spin dynamics in FKLM systems such as manganites and 
DMS, but also allows seamless extention to include correlated band electrons, so that quan- 
tum fluctuation effects associated with both finite S and electron correlation can be studied 
on an equal footing, which is particularly relevant for the ferromagnetic manganites. 
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